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ABSTRACT 

Standing sausage modes in flare loops are important for interpreting quasi- 
periodic pulsations (QPPs) in solar flare lightcurves. We propose an inversion 
scheme that consistently uses their periods P and damping times r to diagnose 
flare loop parameters. We derive a generic dispersion relation governing linear 
sausage waves in pressure-less straight tubes, for which the transverse density 
inhomogeneity takes place in a layer of arbitrary width I and is of arbitrary form. 
We And that P and r depend on the combination of [R/vai, L/R,l/R, pi/pe], 
where R is the loop radius, L is the looplength, vai is the internal Alfven speed, 
and Pi/Pe is the density contrast. For all the density profiles examined, P and 
r experience saturation when L/R S> 1, yielding an inversion curve in the 
[R/vAi,l/R, pi/Pf/\ space with a specific density profile when L/R is sufficiently 
large. When applied to a spatially unresolved QPP event, the scheme yields that 
R/vai is the best constrained, whereas l/R corresponds to the other extreme. 
For spatially resolved QPPs, while L/R S> 1 cannot be assumed beforehand, an 
inversion curve remains possible due to additional geometrical constraints. When 
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a spatially resolved QPP event involves another mode, as is the case for a recent 
event, the fnll set of [fAi, Pi/Pe] can be inferred. We conclnde that the proposed 
scheme provides a usefnl tool for magneto-seismologically exploiting QPPs. 

Subject headings: magnetohydrodynamics (MHD) — Snn: flares — Snn: corona 
— Snn: magnetic helds — waves 


1. INTRODUCTION 

The original ideas that laid the fonndation for the held of solar magneto-seismology 
(SMS) were pnt forward in the 1970s (Uchida 1970, Rosenberg 1970, Zaitsev & Stepanov 1975, 
see also Roberts et ah 1984). However, this held honrished only after a rich variety of low- 
freqnency Magnetohydrodynamic (MHD) waves and oscillations were identihed with the 
advent of the TRACE, SOHO, Hinode, and SDO satellites (for recent reviews, see Banerjee 
et al. 2007, De Moortel & Nakariakov 2012, Mathiondakis et ah 2013; and also Ballester 
et al. 2007, Nakariakov & Erdelyi 2009, Erdelyi & Goossens 2011 for three recent topi¬ 
cal issnes). It is also indispensable to rehne the theoretical nnderstanding of the collective 
wave modes snpported by magnetized tnbes, thereby enabling one to employ the measnred 
wave properties to infer the solar atmospheric parameters that are difhcnlt to measnre di¬ 
rectly (e.g., Roberts 2000; Nakariakov & Verwichte 2005). Regarding its applications to the 
solar corona, SMS can oher snch key information as the magnetic held strength in coronal 
loops (e.g., Nakariakov & Ofman 2001; Erdelyi & Taroyan 2008; Ofman & Wang 2008; White 
& Verwichte 2012) and above streamer stalks (Chen et al. 2010, 2011), the magnitnde of held- 
aligned loop hows (Li et al. 2013; Chen et al. 2014), the temperatnre of loop plasmas (e.g.. 
Marsh & Walsh 2009), the coronal ehective adiabatic index (Van Doorsselaere et al. 2011), 
as well as the longitndinal (Verth & Erdelyi 2008; Andries et al. 2009; Lnna-Cardozo et al. 
2012a) and transverse strnctnring (e.g., Arregni et al. 2007; Goossens et al. 2008; Yuan et al. 
2015). In addition, SMS applications with torsional Alfven waves have proven invaluable 
in inferring the magnetic held structure at chromospheric heights (Jess et al. 2009; Fedun 
et ah 2011). Likewise, Luna-Cardozo et al. (2012a) demonstrated the potential of using 
longitudinal waves to infer the longitudinal variation of density and magnetic held strength 
in chromospheric waveguides. 

Magneto-seismological applications with standing kink modes (with azimuthal wavenum¬ 
ber m = 1) have been a common practice since their detection with TRACE (Aschwanden 
et al. 1999). Kink oscillations tend to experience substantial damping (e.g., Ofman & As¬ 
chwanden 2002; Ruderman & Erdelyi 2009; Verwichte et al. 2013, and references therein), 
which is usually interpreted in terms of resonant absorption (Ruderman & Roberts 2002, 
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Goossens et al. 2002, also Hollweg & Yang 1988, and the comprehensive review by Goossens 
et al. 2011). With this interpretation, Rnderman & Roberts (2002) and Goossens et al. 
(2002) snggested that the measnred period P and damping time r can be used to infer 
the lengthscale I of the density inhomogeneity across coronal loops in units of loop radius 
R. For this purpose, the largely unknown transverse density distribution was shown to be 
important, since its formulation has a considerable impact on P and r (Soler et al. 2013, 

2014) . 

While kink modes have attracted much attention, sausage modes (with m = 0) are 
equally important in SMS. In fact, sausage modes are even more important from the stand¬ 
point of solar atmospheric heating given their stronger compressibility and ubiquity in the 
lower solar atmosphere (Morton et al. 2011, 2012; Freij et al. 2014; Dorotovic et al. 2014; 
Grant et al. 2015; Moreels et al. 2015). In addition, sausage modes are important for inter¬ 
preting quasi-periodic pulsations (QPPs) in the lightcurves of solar flares (see Nakariakov 
& Melnikov 2009, for a recent review). Two distinct regimes are known to exist, depending 
on the axial wavenumber k along flare loops (Nakariakov &: Verwichte 2005). The trapped 
regime results when k exceeds some critical value fee, where the energy of sausage modes 
is well conhned to magnetic tubes. When k < fcc, the leaky regime arises and sausage 
modes experience apparent temporal damping by radiating their energy into the surround¬ 
ing fluid (Spruit 1982; Gaily 1986). It is known that k^ depends sensitively on the density 
contrast between loops and their surroundings (e.g., Kopylova et al. 2007). In addition, both 
eigen-mode analyses (Kopylova et al. 2007; Vasheghani Farahani et al. 2014) and numerical 
simulations from an initial-value-problem perspective (Nakariakov et al. 2012; Ghen et al. 

2015) indicated that the period P of sausage modes increases smoothly with decreasing k (or 
equivalently with increasing looplength L given that k = n/L for fundamental modes) until 
reaching some saturation value Ps for sufficiently thin loops {R/L 1). Likewise, identically 
inhnite in the trapped regime for ideal MHD fluids, the attenuation time r decreases with 
decreasing k before experiencing saturation at Tg when R/L 1. 


Magneto-seismological applications of sausage modes are possible due to their depen¬ 
dence on atmospheric parameters (Luna-Gardozo et al. 2012a,b). The practice based on the 
measured period and damping time can be illustrated by the study presented in Kopylova 
et al. (2007) where a step-function (top-hat) form was adopted for the transverse density 
distribution. The saturation values, Pg and Tg/Pg, for large density contrasts are approxi¬ 
mately 
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As an example, Kopylova et al. (2007) examined the QPP in the radio emissions reported 
in McLean & Sheridan (1973), where P and r/P were found to be ~ 4.3 secs and ~ 10, 
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respectively. With the damping attributed to wave leakage, Eq. (1) then yields a density 
contrast pjPe of ~ 100, and a transverse Alfven transit time R/vm of ~ 1.64 secs, provided 
that the flaring loop in question is sufficiently thin. However, the dependence of P and 
r on loop parameters is substantially more involved if one goes a step closer to reality by 
replacing a step-function density prohle with a smooth one. Even for thin loops where neither 
P nor r depends on looplength, both Pg (Nakariakov et ah 2012) and Tg (Chen et al. 2015) 
may be sensitive to the steepness, or equivalently the lengthscale, of the transverse density 
distribution. In mathematical terms, this means that for a given density prohle. 



( 2 ) 


When only Pg and Tg are known as is the case for measurements without imaging capabilities, 
the appearance of //P no longer allows a unique pair of [pi/pe, P/vAi] to be deduced. Despite 
this, one can still constrain the combination [R/vai, Pi/PeJ/R] by developing a scheme in 
much the same way that kink modes were employed (Arregui et al. 2007; Goossens et al. 
2008; Soler et al. 2014), the only difference being that the transverse Alfven time R/va\ 
replaces the longitudinal one L/vai- 

The present study aims to develop the aforementioned scheme employing measured 
periods and damping times of standing sausage modes. An essential ingredient will be 
to establish the functions P and G in Eq. (2). To this end, we will derive an analytical 
dispersion relation (DR) governing linear sausage waves hosted by magnetized tubes with 
a rather general transverse density distribution. The only requirement here is that this 
density distribution can be decomposed into a uniform cord, a uniform external medium, 
and a transition layer connecting the two. However, the density distribution in the transition 
layer is allowed to be arbitrary, thereby making the DR applicable to a rich variety of density 
prohles. We note that this kind of density prohles has been extensively adopted in kink mode 
studies (e.g., Soler et al. 2014, and references therein). We further note that developing an 
analytical DR is important in its own right. Apart from the step-function prohle (e.g.. Spruit 
1982; Cally 1986), analytical DRs in the cylindrical case are available only for a limited set 
of density prohles (Edwin &: Roberts 1986, 1988; Lopin &: Nagorny 2014). 

This manuscript is organized as follows. Section 2 presents the derivation of the DR and 
our solution method. A parameter study is presented in Sect. 3 to examine how the period 
and damping time of sausage modes depend on tube parameters, thereby establishing our 
numerical scheme for inverting measurements of spatially unresolved QPPs. An extension 
to spatially resolved QPPs is then given in Sect. 4. Finally, Sect. 5 closes this manuscript 
with our summary and some concluding remarks. 
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2. MATHEMATICAL FORMULATION 
2.1. Description for the Equilibrium Tube 

We consider sausage waves in a structured corona modeled by a density-enhanced cylin¬ 
der with radius R aligned with a uniform magnetic held B = Bz, where a cylindrical coor¬ 
dinate system (r, 9, z) is adopted. The equilibrium density is assumed to be a function of r 
only and of the form 


p(r) 


Pi, 0 < r < n = i? — //2, 
Ptr{r), ri<r <re = R + l/2, 

Pe, r > re. 


( 3 ) 


The prohle between [rj, is such that the equilibrium density p decreases continuously from 
the internal value p\ to the external one pe- The thickness of this transition layer, denoted 
by /, is bounded by 0 and 2R. The former represents the steepest prohle of a step-function 
form, whereas the latter corresponds to the least steep case. 

While our analysis is valid for arbitrary prescriptions of ptr, a number of choices have to 
be made to evaluate quantitatively the ehects of equilibrium density prohles. To this end, 
we select the following prohles. 


ptrir) 
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pi ■ 

Pe 
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Pi ~ Pe 
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linear, 
parabolic, 
inverse — parabolic, 
sine. 


( 4 ) 


The prohles labeled linear, parabolic and sine have been examined in substantial detail in 
the context of standing kink modes (Soler et al. 2013, 2014). An additional prohle, labeled 
inverse-parabolic, is added to make the list more comprehensive in that it naturally com¬ 
plements the parabolic one. Figure 1 illustrates the r-dependence of the chosen equilibrium 
density prohles. For illustration purposes, pi/pe is chosen to be 100, and l/R\s chosen to be 
unity. 
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2.2. Solutions for Radial Lagrangian Displacement and Total Pressure 


Perturbation 


Appropriate for the solar corona, we work in the framework of cold (zero-/?) MHD, in 
which case sausage waves do not perturb the 2 ;-component of the plasma velocity. Let 5vr 
denote the radial velocity perturbation, and let and denote the radial and longitudinal 
components of the perturbed magnetic held hb, respectively. The perturbed total pressure 
is then 5ptot = B ■ (5b/47r = BShz/ATT given the absence of thermal pressure in the zero-/? 
limit. Fourier-decomposing any perturbed value 6f{r,z;t) as 



(5) 


one hnds from linearized, ideal, cold MHD equations that 



( 6 ) 


where the prime ' = d/dr. In addition, = ivrjuj is the Fourier amplitude of the radial La¬ 


grangian displacement, and vpSf) = B/^J4:^Tp{r) is the Alfven speed. The Fourier amplitude 
of the perturbed total pressure is 



(7) 


With azimuthal wavenumber m being 0, the equations governing linear sausage waves 
are free of singularities, making our derivation simpler than in kink mode studies where a 
treatment of singularity is necessary (Soler et ah 2013). To be specihc, the solutions to 
Eq. (6) in the transition layer can be expressed as a regular series expansion in x = r — R. 
Let ^tr,i and ^tr ,2 denote two independent solutions. 


OO 


OO 



( 8 ) 


Without loss of generality, one may choose [ao,ai] = [R, 0] and [&o,^i] = [0,1]- Expanding 
the equilibrium density p about x = 0 as well, one hnds that 


OO 



(9) 


with po = p\x=o and 


1 R^p{x) 


, n > 1. 


Pn 


n\ dx” 


x=0 


( 10 ) 
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Plugging Eq. (8) into Eq. (6) with the change of independent variable from r to x, and then 
ordering the terms according to powers in x, one hnds that 

X 2 = ^ [RXi + {r]u‘^R^Po - k'^R^ - 1) Xo] , 

X3 = ^ [QRX2 + {vuj^R^Po - k'^R^) Xi 

+ {rjoj'^R^Pi + 2rioj’^RpQ - 2k‘^R) Xo\ , (11) 

-1 f 

- 7^{{n-l){2n-?,)Rxn-i + pi^‘^R^'y^pn-i-2Xi 

n[n — IjR-^ I 

n —3 

+ [{n - 3) (n - 1) - k'^R'^] Xn -2 + 2pu‘^R pn-i-sXi 

1=0 

n —4 ^ 

-2k‘^RXn-3 - k‘^Xn-4 + ^ pn-i-iXi > , {n > 4:) 

1=0 ) 

where p = Air /and x represents either a or b. 

With Eq. (8) at hand, the solution to Eq. (6) can be expressed as 


r AJi(/iir), 0 < r < ri, 

ir{r) = I Aiftr,l(a;) + Mitrpix), T, < T < Te, 
[ r->re, 


( 12 ) 


where Ai and A 2 are arbitrary constants, and Jn and are the n-th-order Bessel 

and Hankel functions of the hrst kind, respectively (here n = 1). In addition, /if^ = 
“ k"^ n^j g = B"^/{Airpi^^). As discussed in Cally (1986), requiring that —7r/2 < 
arg/ii, arg/ie <t^/2, does not exclude any additional independent solution. Furthermore, ex¬ 
pressing the external solntion in terms of Hn'^ permits a unihed examination of both trapped 
and leaky waves. Indeed, the trapped regime arises when arg/ig = vr/2, in which case one 
hnds that H‘^\ppr) = —{2/ir)Ki{ar) with a = p^ji being real and positive (see discussions 
on page 281 in Cally 1986). Now the Fonrier amplitude for the total pressure perturbation 
can be evaluated with Eq. (7), the results being 


Vtot{r) = < 


47r 

52 

Airr 


P{Jq{Pit') , 

d 




X + R)^tr,l{x) 


+ ^2^ 
da; 


[X + R)itr,2{x) 


rrWr I 

-/ieh/o APer) , 


Air 


0 < r < Tj, 

ri < r < re (13) 

r > re. 














The expression for ptot for the ranges r < and r > Ve can also be independently verihed 
by using an alternative relation between and ptot, 


dptot/dr 
p{u'^ — 


2.3. Dispersion Relation of Sausage Waves 


The dispersion relation (DR) governing linear sausage waves can be derived by requiring 
that both and ptot be continuous at the interfaces r = n and r = r^. This leads to 


AhiD'^o(AiTi) 

AePereH^^\pere) 


^l^tr,l(^e) T ^2^tr,2(^e)) 

+ ^21^,2 (a^i)] + + ^2'f(r,2 (^i)] > 

[Aiftr,l(a;e) + ^2|tr,2(Xe)] + Tg ^ (Xg) + ^21^,2(^^e)], 


where Xi = —1/2 and = 1/2. Eliminating Ai (Ag) by dividing the third (fourth) by the 
hrst (second) equation, one hnds that 


AiAi + A 2 A 2 — 0, 

A 3 R 1 + A 4 R 2 = 0 , 


where the coefficients read 


Al 'Ctr,l(^i) T (^i) A^i'Ctr,! (^^i)) 

A 2 = {tr,2(a;i) + riCtrpiXi) - Xiitpixi), 

A3 Ctr,1(2^0) T ^e'btr,l(^e) ^eCtr,l(^e)) 

A4 'Ctr,2(^e) T ^e'btr,2(^e) ^e'Ctr,2(^e)) 


(16) 


with 


Xi 

Xg 


k^iriMpiTj) 

Ji(piri) 

PereHQ\pere) 

H[^\pere) 


(16) 


Evidently, for Eq. (14) to allow non-trivial solutions of [^ 1 ,^ 2 ], one needs to require that 


A 1 A 4 — A 2 A 3 — 0 , 


( 17 ) 


which is the DR governing sausage waves in nonuniform loops. 
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Does Eq. (17), the DR valid for arbitrary l/R, recover the well-known result for the 
step-function prohle when l/R —)■ 0? Retaining only terms to the 0-th order in l/R and 
noting that Ui Ue ~ R, one hnds that A„ (n = 1, • • • , 4) simplify to 

— (1 — 4“ Rci±y 

A 2 = (1 — Xijbo + Rbi, 

A3 = (1 — A'e)ao -l- Rai, 

A 4 = (1 — Xejbo + Rbi- 


Substituting these expressions into Eq. (17), one hnds that 


(Xi — Xe)(ai6o — Oq&i) — 0. 


This leads to X, = Xg, given that aibo — ao&i is not allowed to be zero for and ^tr ,2 to 
be independent. In other words (see Eq. (16)) 

RiJo{RiR) _ {fleR) n 8i 

JlifiiR) ~ H[^\fi^R) 

which is the DR for equilibrium density prohles of a step-function form (e.g.. Gaily 1986). 


While the DR is equally applicable to propagating waves, we will focus on standing 
modes for which the axial wavenumber k is real, while the angular frequency u is allowed to 
be complex-valued. Furthermore, let ns focus on fundamental standing modes supported by 
magnetized tubes of length L. In this case, another measure to validate the DR, independent 
of the eigen-value problem approach, is to employ the linearized, time-dependent, cold MHD 
equations to derive an equation governing the transverse velocity perturbation Vr{r, z, t). In 
view of the line-tying boundary conditions at the loop ends z = 0 and z = L, one may 
express Vr{r,z,t) as v{r,t) sin{kz) with k = tt/L, yielding (Nakariakov et al. 2012; Chen 
et al. 2015) 


d'^v{r,t) 2 / 1 ^ 

df2 A ^ 



(19) 


When supplemented with appropriate boundary and initial conditions, the signal of Vr{r,t) 
at some arbitrarily chosen distance from the tube axis can be followed. As demonstrated 
in Nakariakov et al. (2012); Chen et al. (2015), after a transitory phase this signal evolves 
into a harmonic (decaying harmonic) form when k is larger (less) than some critical value, 
corresponding to the well-known trapped (leaky) regime. Numerically htting the signal 
with a sinusoidal (exponentially decaying sinusoidal) function then yields the period P {P 
together with the damping time r) for trapped (leaky) modes. We also adopt this approach 
and compare the derived values for P and r with what is found by solving the DR for 
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complex-valued cj at given real k. As will be shown in Fig. 2, the two sets of independently 
derived values for [P, r] agree remarkably well. At this point, it should be remarked that 
once a choice for ptr is made, P and r depend only on the combination of parameters 
[pi/pe, l/R, L/R] when they are measured in units of the internal Alfven transit time R/va\- 
Here we have used L/R = 7r/{kR) in place of the dimensionless axial wavenumber kR. 

In general, the DR (Eq. (17)) is not analytically tractable and is solved numerically 
for a given ptr prohle and some given combination of [pi/Pe,l/R, L/R], To do so requires 
the inhnite series expansion in Eq.(8) to be truncated by retaining the terms with n up to a 
certain N. A value of A^ = 101 is chosen for all the numerical results to be presented, and we 
have made sure that using an even larger N does not introduce any appreciable difference. 
In addition, we focus only on the lowest order modes, namely those with the simplest radial 
structure in the eigen-functions. When verifying these results with the computations from 
an initial-value-problem perspective, we choose not to use a too localized initial perturbation 
v{r,t = 0) (see Eq.(19)), otherwise higher order modes are introduced to contaminate the 
v{r,t) signals. 


3. NUMERICAL RESULTS AND THE INVERSION SCHEME 

Let us start with an examination of how the ptr prohle impacts the dispersive properties 
of standing sausage modes. Figure 2 shows the dependence on the length-to-radius ratio 
L/R oi the period P and damping time r for different choices of ptrif) as labeled. For 
illustration purposes, we choose the density contrast p-JPe to be 100, and choose l/R to 
be unity. The black line in Fig. 2a, which represents 2L/v\e, separates the trapped (to 
its left, where r is identically inhnite) from leaky (to its right) regimes. The curves are 
found by solving the analytical DR (Eq. (17)), whereas the open circles are found by solving 
the corresponding time-dependent equation (see Eq. (19) and the associated description). 
Evidently, the periods P and damping times r obtained from the two independent approaches 
agree with each other remarkably well. Figure 2 indicates that the overall tendency for P (r) 
to increase (decrease) with L/R is seen for all the equilibrium density prohles considered. 
In particular, regardless of the prohles, both P and r tend to some asymptotic values at 
large L/R. However, the choice of equilibrium density prohles has a considerable inhuence 
on the specihc values for P and r. This is particularly true if one compares the results for 
the parabolic and inverse-parabolic prohles, given by the green and blue curves, respectively. 
Furthermore, while the periods P for the linear and sine prohles diher little for the chosen 
l/R, the damping times r show a stronger prohle dependence. This signihes the importance 
of using P and r in a consistent manner when one attempts to deduce how the transverse 
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equilibrium density is structured. 

The effects of equilibrium density profile can be better brought out by capitalizing on 
the fact that both P and r experience saturation for sufficiently thin loops. Let and 
Ts denote the saturation values. Figure 3 presents how Pg and Tg/Pg depend on //P, the 
width of the transition layer in units of loop radius. The density contrast is chosen to be 
100, and different choices of ptr are represented by the curves in different colors as labeled 
in Fig. 3b. In addition, in place of Tg, the ratio Tg/Pg is plotted since it is a better measure 
of signal quality. One sees from Figs. 3a and 3b that the curves converge at //P —)■ 0 as 
expected given that the DR (Eq. (17)) simplifies to Eq. (18) pertinent to a top-hat density 
distribution. Figure 3a indicates that the //P-dependence of Pg critically depends on how 
Ptr is described. For the parabolic profile, Pg increases monotonically with I /P, whereas the 
opposite trend is found for the inverse-parabolic profile. When it comes to the linear and 
sine profiles. Fig. 3a shows that the //P-dependence of Pg is not as strong, and the difference 
between the two profiles is discernible only when //P > 1. Moving on to Fig. 3b, one sees 
that the ratio Tg/Pg decreases monotonically with l/R for all the profiles, meaning that wave 
leakage plays an increasingly important role in attenuating sausage modes when the loop 
becomes more diffuse. Reinforcing the impression from Fig. 2b, one sees that relative to Pg, 
Tg/Pg better discriminates the equilibrium density profiles. 

So far we have fixed the density contrast p-JPe at 100. One naturally asks what happens 
if Pi/Pe is varied? Figure 4 presents the distribution of Pg (the left column) and Tg/Pg (right) 
in the [pi/pe,l/R] plane. Each row represents one of the four density profiles as labeled. 
Besides, the red curve represents where Tg/Pg = 10, the value for the QPP event reported 
in McLean &: Sheridan (1973). Examine the left column first. One sees from Figs. 4b and 
4c that in the parameter range examined, Pg tends to increase (decrease) with l/R aX any 
given Pi/pe when the parabolic (inverse-parabolic) profile is chosen. In contrast. Figs. 4a and 
4d indicate that Pg for the linear and sine profiles shows a nonmonotonical dependence on 
l/R, even though this variation is hardly discernible. Now consider the right column, from 
which one can see that regardless of the profiles, Tg/Pg decreases with increasing l/R for all 
the density contrasts examined. In addition, the dependence of Tg/Pg on //P is the strongest 
for the inverse-parabolic profile, and the least strong for the parabolic one. The dependence 
for the linear and sine profiles lies in between, with the dependence in the linear case being 
slightly stronger. 

Conceptually, Fig. 4 can be used to invert the measured values of the period and damping 
time of sausage modes, provided that the loops hosting these oscillations are sufficiently thin. 
Consider the QPP event reported in McLean & Sheridan (1973) as an example, for which 
Pg = 4.3 secs and Tg/Pg = 10. It then follows that for a given density profile, any point 
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along the corresponding red curve in Fig. 4 can equally reproduce the measured Tg/Ps- After 
reading any pair of [pi/pe,l/R], one can read from the left column the corresponding value 
for Pg in units of R/va\. With Pg known, one can then deduce R/va\. In practice, however, 
constructing such a contour plot is not necessary, and one may simply consider the following 
3-step inversion scheme. First, one starts with the dispersion relation for a step-function 
density profile (//i? = 0), Eq. (18), to find the value for the density contrast p-JPe such that 
Tg/Pg agrees with the measured value. Second, with the p\/Pe value for a smaller l/R as a. 
good guess, one can then solve Eq. (17) to hnd a new pi/Pe that yields the measured Tg/Pg 
by increasing l/R from 0 to 2 consecutively. Third, with p-Jp^ found for all possible l/R, one 
can solve Eq. (17) for a given pair of [pi/Pe,l/R], yielding a value for Pg/{R/vxi). Finding 
the transverse Alfven transit time R/v^i is then straightforward since Pg is known. 

The product of the inversion scheme is an inversion curve in the three-dimensional (3D) 
space formed by R/va\, l/R and pi/ Pe- Figure 5 presents such curves (the solid lines) and 
their projections onto various planes (dashed) for the examined density profiles, pertinent 
to the event reported in McLean &: Sheridan (1973). To help digest this hgure, a number 
of points are read from the curves and presented in Table 1. One sees that among the 
parameters forming this 3D space, R/vai is the best constrained. The biggest (smallest) 
value, 2.13 (1.18) secs, is found for the inverse-parabolic (parabolic) profile when l/R ^ 2. 
In other words, the biggest value exceeds the smallest one by only 79.7%. As to the density 
contrast, the biggest value (251.4 found for the inverse-parabolic profile when l/R ^ 2) is 
larger than the smallest one (88.1 when //i? —)■ 0) by 185%. The least constrained parameter 
is l/R, with any value in the allowed range from 0 to 2 being possible. 


4. FURTHER DEVELOPMENT OF THE INVERSION SCHEME 

Before proceeding, let us first recap the key points in the scheme for inverting the 
measured period P and damping time r of sausage modes. From the outset, we have assumed 
that only P and r are known, whereas we have no information on either the geometric 
parameters {R, L, and 1) or the physical parameters (uai, Ph pg). On top of that, 
the specihc density profile is also assumed to be unknown. This happens when one has only 
spatially unresolved observations (see e.g., the majority of the events compiled in Aschwanden 
et al. 2004, Table 1). In this case, the dispersion relation, Eq. (17), suggests that in general 
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the periods and damping times can be formally expressed as 



saus 


( 20 ) 


With only two measured values available, any point on a two-dimensional (2D) surface 
in the four-dimensional (4D) space formed by [R/vpa, L/R,l/R, p-J can reproduce the 
measurements, even if one is allowed to prescribe a density prohle. If the sausage modes are 
in the trapped regime in the sense that the QPPs do not show temporal damping, then the 
situation becomes even less desired since now the restriction from r is no longer available. 
This complexity can be alleviated if the condition L/R ^ 1 holds for the flare loops in 
question since L/R \io longer appears such that Eq. (2) is restored. One then Ends a curve 
in the [R/vm^I/R, p-Jpe] space as shown by Fig. 5. 

The situation improves if the QPP events are spatially resolved, since the looplength 
(L) and the outer interface of the loops (re = R + 1/2) can be considered known. In this 
case L/R and l/R are no longer independent but are related by 



( 21 ) 


where we have used the relation 


( 22 ) 


2 +l/R' 


Now the situation is similar to what Eq. (2) implies: if the sausage mode is a trapped (leaky) 


one, then a 2D surface (ID curve) can be deduced in the [R/vpj^^^l/R, p-Jp^ space. It is just 
that now for a point on this 2D surface or along the ID curve, one can further deduce ua) in 
view of Eq. (22). 

Something fascinating happens if one observes a spatially resolved QPP event hosting 
more than just a sausage mode. For illustration purposes, we consider the situation where a 
temporally decaying kink mode is involved, and its damping can be attributed to resonant 
absorption. In this case, the period and damping time for the kink mode can be formally 
expressed by 



( 23 ) 
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Note that while L/R m principle can be incorporated into Eq. (23), in reality there is no 
need to do so since the corrections to Pkink and Tkink due to hnite L/R are of the order 
{R/LY (Van Doorsselaere et ah 2004; Goossens et ah 2008). Even for relatively thick flare 
loops, R/L is of the order 0.1 and these corrections amount to only a few percent. Note 
further that Pkink and Gkink have been extensively studied, and a graphical representation 
can be found in Fig. 1 of Soler et ah (2014). With Pkink, Gkink, Psaus, and Gsaus known, it 
is then possible to fully constrain the unknowns [/, P, UAi, Pi/Pe], if one assumes a density 
prohle. In fact, the measured values of [Pkink, PLink, Psaus, Pans] are more than sufficient: 
I and R are not independent but are related through Eq. (22). This suggests that only 
three expressions contained in Eqs. (20) and (23) are needed. In practice, we consider the 
expression for the kink mode period as the redundant one. Now the inversion procedure is 
rather straightforward. In view of the relation (21), the two expressions for the damping- 
time-to-period ratio in Eqs. (20) and Eqs. (23) contain only two unknowns, namely l/R 
and Pi/pe. With both Tsaus/Psaus and Tkink/Pkink available, one can deduce a unique pair 
of [//P, pi/pe]- The loop radius P then follows from Eq. (22), which then enables one to 
evaluate va\ with the first expression in Eq. (20). Finally, as a safety check, one can proceed 
to evaluate, with the first expression in Eq. (23), the theoretically expected kink mode period 
Pkink,theory The deviation of Pkink,theory from the measured one then allows to say a few words 
on how safe it is to identify the oscillation signals with some particular modes. 

As an illustration of the aforementioned inversion procedure, let us consider the QPP 
event in microwave emissions measured with the Nobeyama Radioheliograph (NoRH) on 
14 May 2013 (Kolotkov et ah 2015). Lucky enough, it is likely that this event contains 
a fundamental kink mode with Pkink = 100 secs and Tkink/Pkink = 2.5, in addition to a 
fundamental sausage mode with Pgaus = 15 secs and Tgaus/Psaus = 6. Assuming that the 
apparent width measured therein corresponds to 2re, one hnds that Te = 4 x 10^ km and 
L = 4 X 10^ km, meaning that Ljr^ = 10. Let us assume that wave leakage is responsible 
for damping the sausage mode, and resonant absorption is responsible for damping the kink 
one. Furthermore, let us assume that the sine profile best describes the equilibrium density 
distribution. In this case, the analytical expressions obtained in the thin-tube-thin-boundary 
approximation are accurate to within ~ 25% (e.g.. Van Doorsselaere et al. 2004; Soler et ah 
2014). This enables us to illustrate our inversion procedure without resorting to a fully 
numerical solver to establish Pkink and Gkink- Now the formal expressions given by Eq. (23) 
can be replaced with (e.g., Goossens et al. 2008) 



T / 2(1 + Pi/pe) 
^AiV P/Pe 


Tkink 


2 pi/Pe + l 1 

7rpi/pe - 1 


( 24 ) 
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Following the outlined inversion procedure, we find that l/R = 0.272 and pi/pe = 29.8 as 
constrained by the ratios r/P. It then follows that R = 3.52 x 10^ km, from which one can 
deduce that v^i = 623 km s“^. Finally, Fkink,theory is found to be 92.2 secs, which agrees 
with the measured value of 100 secs to within 8%. This safety check lends support to the 
interpretation of the two modes in terms of fundamental kink and sausage modes as done 
by Kolotkov et ah (2015). Besides, the deduced Alfven speed vm and density contrast p-JPe 
both seem reasonable. On top of that, there is no need for one to worry too much about 
the accuracy of Eq. (24) describing the kink mode: at the deduced [//i?, pi/pe], this equation 
yields P and t/P to an accuracy better than 1% and ~ 6%, respectively (see Fig. 1 in Soler 
et al. 2014). 


5. CONCLUSIONS 

How plasma density is structured across various magnetic structures in the solar corona 
remains largely unknown. It has been a common practice to deduce this key information by 
employing magneto-seismological techniques that invert the measured period P and damping 
time r of standing kink modes collectively supported by a magnetic structure (e.g., Goossens 
et al. 2008; Soler et al. 2014). In contrast, while quasi-periodic pulsations (QPPs) in the 
lightcurves of solar flares are often attributed to standing sausage modes in flare loops and 
therefore can also offer the associated period and damping time, a scheme is missing for 
inverting these two measurements to deduce the information on the density distribution 
transverse to flare loops. The primary aim of this study has been to construct such a 
scheme. To this end, we worked in the framework of cold (zero-/3) MHD and modeled flare 
loops as straight cylinders with a transverse density profile characterized by a transition 
layer sandwiched between a uniform cord and a uniform external medium. An analytical 
dispersion relation (DR) governing linear sausage waves, Eq. (17), was derived by solving the 
perturbation equations in terms of a regular series expansion in the transition layer. This 
DR, valid for arbitrary choices of the density profile in the transition layer, then enabled us 
to examine the effects of density structuring on the periods and damping times of sausage 
modes, thereby facilitating the construction of the inversion scheme. 

In general, we found that P and r of sausage modes depend on a combination of param¬ 
eters [P/vAi, L/R,l/R, Pi/pe] as formally expressed by Eq. (20), where the functions Fsaus 
and Gsaus are a product of the DR. Here L and R denote the looplength and loop radius, 
respectively. Furthermore, I is the width of the transition layer, uai is the Alfven speed in 
the cord, and pi/pe is the density contrast between the loop and its surrounding fluid. We 
showed that for the density profiles examined, both P and r experience saturation for suffi- 
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ciently large L/R when the rest of the four parameters are hxed. The choice of the transverse 
density prohle was found to have a considerable influence on P and r, their dependence on 
l/R in particular. 

Our inversion scheme can hnd applications to both spatially unresolved and resolved 
QPP events. For spatially unresolved ones, we showed that the best one can do is to deduce 
a ID curve in the [R/vAi,l/R, pi/pe] space for a prescribed density prohle. This happens if 
the QPPs in question experience temporal damping, and if the hare loops hosting them can 
be assumed to be sufficiently thin. When applied to a QPP event reported by McLean & 
Sheridan (1973), this inversion technique indicated that the transverse Alfven transit time 
R/vAi is the best constrained, varying by a factor of 80% if the uncertainties in specifying the 
density prohle are taken into account. The density contrast Pi/Pe is less well constrained, 
with the largest deduced value exceeding the smallest one by a factor of 1.85. The least 
constrained is the transverse density lengthscale in units of loop radius {l/R), any value in 
the allowed range (0, 2) can be equally possible to reproduce the measurements. 

When it comes to spatially resolved events, the geometric parameters L and R + l/2 are 
additional constraints, on top of the measured values for P and r. Even though in this case 
one cannot assume L/R ^ 1 a priori, it is possible to deduce a ID curve in the [ uai , I, Pi/Pe] 
space for a chosen density prohle since R is expressible in terms oi l/R. If a spatially resolved 
QPP event comprises more than just one sausage mode, then it is possible to deduce the 
full information on [I, R,VAi, Pi/Pe]- In this case, the inversion problem becomes an over¬ 
determined one. We illustrated this fascinating application with the QPP event reported 
by Kolotkov et al. (2015) where a fundamental kink and a fundamental sausage mode were 
suggested to co-exist and both experience temporal damping. Attributing the temporal 
damping of the kink mode to resonant absorption, and that of the sausage mode to wave 
leakage, we were able to deduce the full set of [I, R,VAi, Pi/Pe] by using Eqs. (20) and (24) 
and assuming a sine prohle for the density distribution. One redundant equation, here taken 
as the expression for the kink mode period, can then allow a safety check on, say, whether 
it is reasonable to interpret the signals in the QPP event as the aforementioned modes. For 
this particular event, our results demonstrated that the interpretation provided in Kolotkov 
et al. (2015) is reasonable, not only because the deduced parameters seem realistic, but also 
because the theoretical prediction for the kink mode period agrees with the measured one 
remarkably well. 

Our scheme nonetheless has a number of limitations. First, we have assumed that 
the temporal damping of sausage modes is due to wave leakage, an ideal MHD process. 
In reality, non-ideal mechanisms like electron heat conduction and ion viscosity can provide 
additional channels for damping sausage modes. While these non-ideal processes were shown 
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by Kopylova et al. (2007) to be unlikely the cause for the temporal damping in the QPP 
event reported by McLean & Sheridan (1973), their importance has yet to be assessed for 
the event reported by Kolotkov et al. (2015). Second, working in the cold MHD, we have not 
taken into account the possible effects due to hnite plasma beta, which may be of the order 
unity in flare loops. However, the corrections due to hnite beta seem marginal (Inglis et al. 
2009). Third, the longitudinal variation in neither the plasma density nor the magnetic held 
strength has been considered, even though the corrections due to this variation are unlikely 
to be signihcant (Pascoe et al. 2009). Fourth, the density inhomogeneity in hare loops was 
assumed to be in a monolithic form, whereas in reality these loops may be multi-stranded. 
While the hne structuring in the form of randomly distributed concentric shells is found to 
have a far less signihcant inhuence than the monolithic component of the density distribution 
(Chen et al. 2015, see also Pascoe et al. 2007), there is a need to rigorously assess the ehects 
due to hne structuring in the form of randomly distributed strands. 


We thank the referee for the constrnctive comments, which helped improve this manuscript 
substantially. This work is supported by the National Natnral Science Foundation of China 
(41174154, 41274176, and 41474149), and by the Provincial Natural Science Foundation of 
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r/R 

Fig. 1.— Transverse equilibrium density profiles as a function of r. The profiles differ only 
in how they are described in a transition layer sandwiched between the internal (with a 
uniform density p\) and external (with a uniform density pe) parts. The transition layer is 
of width /, and is located between r\ = R — l/2 and = R + 1/2, with R being the cylinder 
radius. Four different choices of the density prohles in the transition layer are adopted as 
labeled, and are given by Eq. (4). For illustration purposes, I is chosen to be R, and Pi/pe is 
chosen to be 100. 
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Fig. 2.— (a) Period P and (b) damping time r as functions of loop length L. Four 

different choices of density prohles are examined as labeled. The black line in (a) represents 
P = 2L/vAe and separates the trapped (to its left) from leaky (right) regimes. The open 
circles represent the values for P and r obtained by solving Eq. (19) from an initial-value- 
problem perspective, which is independent from the eigen-value-problem approach presented 
in the text. Here the width of the transition layer / = R, and the density contrast p\/pe = 100. 
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Fig. 3.— Saturation values for (a) period Pg and (b) damping time Tg as functions of 
the width of the transition layer 1. These saturation values are attained for sufficiently thin 
loops. Four different choices of the density prohles are examined as labeled. Here the density 
contrast p-JPe = 100. 
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Fig. 4.— Contour plots in the [1/R, p-J space for the saturation values for period (the 
left column) and damping-time-to-period ratio Tg/Ps (right). These saturation values are 
attained for sufficiently thin loops. Each row represents one of the four different choices of 
the density prohles as labeled. The red curve in each panel represents where Tg/Pg = 10, 
corresponding to the event reported in McLean & Sheridan (1973) 
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Fig. 5.— Inversion curves (the solid lines) together with their projections (dashed) in the 
three-dimensional space formed by [R/vAhl/R, Pi/P<^- Four choices for the density prohles 
are examined and given in different panels as labeled. All points along an inversion curve are 
equally compatible with the quasi-periodic-pulsation event reported in McLean & Sheridan 
(1973) where the oscillation period is 4.3 secs, and the damping-time-to-period ratio is 10. 
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l/R 

linear 

parabolic 

inverse- 

parabolic 

sine 

Pi/Pe 

R/vAi (sec) 

Pi/Pe R/vAi (sec) 

Pi/pe R/vAi (sec) 

Pi/pe R/va\ (sec) 

0.01 

88.1 

1.62 

88.2 

1.62 

88.1 

1.62 

88.1 

1.62 

0.2 

88.9 

1.62 

89.2 

1.57 

89.3 

1.68 

89.1 

1.62 

0.4 

94.6 

1.63 

92 

1.52 

93.2 

1.74 

91.9 

1.62 

0.6 

102.4 

1.63 

95.8 

1.48 

100.4 

1.81 

96.4 

1.63 

0.8 

112.75 

1.62 

100.4 

1.43 

111.5 

1.88 

102.4 

1.63 

1.0 

124.7 

1.61 

105.4 

1.39 

127.5 

1.95 

109.8 

1.62 

1.2 

137.4 

1.6 

110.5 

1.35 

148.4 

2.02 

118.2 

1.62 

1.4 

149.8 

1.57 

115.4 

1.3 

173.2 

2.07 

127.4 

1.61 

1.6 

161.5 

1.54 

120 

1.26 

200.1 

2.1 

136.9 

1.6 

1.8 

172 

1.51 

124.4 

1.22 

226.8 

2.12 

146.4 

1.58 

1.99 

181.1 

1.47 

128.2 

1.18 

251.4 

2.13 

155.4 

1.56 


Table 1: The inverted values for the transverse density length scale in units of loop radius 
(l/R), density contrast pi/pe and the transverse Alfven transit time R/vai- This inversion is 
made for the quasi-periodic-pulsation event reported in McLean & Sheridan (1973), assuming 
that the associated flare loop is sufficiently thin. For this event, the oscillation period is 4.3 
secs, and the damping-time-to-period ratio is 10. 











